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Abstract 

Capacity scaling laws are analyzed in an underwater acoustic network with n regularly located nodes on a 
square. A narrow-band model is assumed where the carrier frequency is allowed to scale as a function of n. In the 
network, we characterize an attenuation parameter that depends on the frequency scaling as well as the transmission 
distance. A cut-set upper bound on the throughput scaling is then derived in extended networks. Our result indicates 
that the upper bound is inversely proportional to the attenuation parameter, thus resulting in a highly power-limited 
network. Interestingly, it is seen that unlike the case of wireless radio networks, our upper bound is intrinsically 
related to the attenuation parameter but not the spreading factor. Furthermore, we describe an achievable scheme 
based on the simple nearest-neighbor multi-hop (MH) transmission. It is shown under extended networks that the 
MH scheme is order-optimal as the attenuation parameter scales exponentially with y/n (or faster). Finally, these 
scaling results are extended to a random network realization. 
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I. Introduction 

A pioneering work of [1], introduced by Gupta and Kumar, characterized the sum throughput scaling 
in a large wireless radio network. They showed that the total throughput scales as 0(-^/n/logn) when 
a multi-hop (MH) routing strategy is used for n source-destination (S-D) pairs randomly distributed in 
a unit arealll MH schemes are then further developed and analyzed in [3]-[9], while their throughput 
per S-D pair scales far slower than 0(1). Recent results [10], [11] have shown that an almost linear 
throughput in the radio network, i.e. 6(n 1_e ) for an arbitrarily small e > 0, which is the best we can 
hope for, is achievable by using a hierarchical cooperation strategy. Besides the schemes in [10], [11], 
there have been other studies to improve the throughput of wireless radio networks up to a linear scaling 
in a variety of network scenarios by using novel techniques such as networks with node mobility [12], 
interference alignment [13], and infrastructure support [14]. 

Along with the studies in terrestrial radio networks, the interest in study of underwater networks 
has been growing with recent advances in acoustic communication technology [15]— [18]. In underwater 
acoustic communication systems, both bandwidth and power are severely limited due to the exponential 
(rather than polynomial) path-loss attenuation with propagation distance and even frequency-dependent 
attenuation. This is a main feature that distinguishes underwater systems from wireless radio links. Hence, 
the system throughput is affected by not only the transmission distance but also the useful bandwidth. 
Based on these characteristics, network coding schemes [17], [19], [20] have been presented for underwater 
acoustic channels, while network coding showed better performance than MH routing in terms of reducing 
transmission power. MH networking has further been investigated in other simple but realistic network 
conditions that take into account the practical issues of coding and delay [21], [22]. 

One natural question is what are the fundamental capabilities of underwater networks in supporting a 
multiplicity of nodes that wish to communicate concurrently with each other, i.e., multiple S-D pairs, 
over an acoustic channel. To answer this question, the throughput scaling for underwater networks was 
first studied [23], where n nodes were arbitrarily located in a planar disk of unit area [1] and the carrier 
frequency was set to a constant independent of n. That work showed an upper bound on the throughput of 
each node based on the physical model assumption in [1]. This upper bound scales as n -V a e -WoiQ{n a )) ^ 
where a corresponds to the spreading factor of the underwater channel, and W represents the branch zero 
of the Lambert function [24] o Since the spreading factor typically has values in the range 1 < a < 2 [23], 
the throughput per node decreases almost as 0(n~ l l a ) for large enough n, which is considerably faster 
than the Q(y/n) scaling characterized for wireless radio settings [1]. 

In this paper, a capacity scaling law for underwater networks is analyzed in extended networks [4], [5], 
[10], [25], [26] of unit node density. Especially, we are interested in the case where the carrier frequency 
scales as a certain function of n in a narrow-band model. Such an assumption changes the scaling behavior 
significantly due to the attenuation characteristics. We aim to study both an information-theoretic upper 
bound and achievable rate scaling while allowing the frequency scaling with n. 

We explicitly characterize an attenuation parameter that depends on the transmission distance and also 
on the carrier frequency, and then identify fundamental operating regimes depending on the parameter. 
For networks with n regularly distributed nodes, we derive an upper bound on the total throughput 
scaling using the cut-set bound. Our upper bound is based on the characteristics of power-limited regimes 
shown in [10], [27]. In extended networks, it is shown that the upper bound is inversely proportional 
to the attenuation parameter. This leads to a highly power-limited network for all the operating regimes, 
where power consumption is important in determining performance. Interestingly, it is seen that unlike 
the case of wireless radio networks, our upper bound heavily depends on the attenuation parameter but 

'We use the following notations: i) fix) = 0(g(x)) means that there exist constants C and c such that fix) < Cg(x) for all x > c. 
ii) f(x) = o(g(x)) means that lim {g = 0. iii) f(x) = n(g(x)) if g(x) = 0(f(x)). iv) f{x) = u(g(x)) if g{x) = o{f{x)). v) 
f(x) = Q(g(x)) if f{x) = 0(g(x)) and g(x) = 0(f(x)) [2]. 

2 The Lambert W function is defined to be the inverse of the function z = W(z)e w ^ z ' and the branch satisfying W(z) > —1 is denoted 

byW (z). 
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not on the spreading factor (corresponding to the path-loss exponent in wireless networks). In addition, to 
constructively show our achievability result for extended regular networks, we describe the conventional 
nearest-neighbor MH transmission [1] with a slight modification, and analyze its achievable throughput. 
It is shown under extended networks that the achievable rate based on the MH routing scheme matches 
the upper bound within a factor of n with arbitrarily small exponent as long as the attenuation parameter 
increases exponentially with respect to *Jn (or faster). Furthermore, a random network scenario is also 
discussed in this work. It is shown under extended random networks that the conventional MH-based 
achievable scheme is not order-optimal for any operating regimes. 

The rest of this paper is organized as follows. Section [TT] describes our system and channel models. In 
Section HHl the cut-set upper bound on the throughput is derived. In Section [TV] achievable throughput 
scaling is analyzed. These results are extended to the random network case in Section |V] Finally, Section 
|VH summarizes the paper with some concluding remarks. 

Throughout this paper the superscript H, f-]^, and || • || 2 denote the conjugate transpose, the [k,i)-\h 
element, and the largest singular value, respectively, of a matrix. I n is the identity matrix of size n xn, tr(-) 
is the trace, and det(-) is the determinant. C is the field of complex numbers and E[] is the expectation. 
Unless otherwise stated, all logarithms are assumed to be to the base 2. 

II. System and Channel Models 

We consider a two-dimensional underwater network that consists of n nodes on a square with unit 
node density such that two neighboring nodes are 1 unit of distance apart from each other in an extended 
network^!, i.e., a regular network [25], [26]. We randomly pick a matching of S-D pairs, so that each 
node is the destination of exactly one source. We assume frequency-flat channel of bandwidth W Hz 
around carrier frequency /, which satisfies / 3> W, i.e., narrow-band model. This is a highly simplified 
model, but nonetheless one that suffices to demonstrate the fundamental mechanisms that govern capacity 
scaling. Assuming that all the nodes have perfectly directional transmissions, we also disregard multipath 
propagation, and simply focus on a line-of-sight channel between each pair of nodes used in [10], [11], 
[27]. Each node has an average transmit power constraint P (constant), and we assume that the channel 
state information (CSI) is available at all receivers, but not at the transmitters. It is assumed that each 
node transmits at a rate T(n)/n, where T{n) denotes the total throughput of the network. 

Now let us turn to channel modeling. An underwater acoustic channel is characterized by an attenuation 
that depends on both the distance r ki between nodes i and k (i,k G {1, • • • ,n}) and the signal frequency 
/, and is given by 

A(r ki ,f) = c r^a(fY k ' (D 

for some constant c > independent of n, where a is the spreading factor and a(f) > 1 is the 
absorption coefficient [16]. For analytical convenience, we assume that the spreading factor a does not 
change throughout the network, i.e., that it is the same from short to long range transmissions, as in 
wireless radio networks [1], [4], [10]. The spreading factor describes the geometry of propagation and is 
typically 1 < a < 2 — its commonly used values are a = 1 for cylindrical spreading, a = 2 for spherical 
spreading, and a = 1.5 for the so-called practical spreading. Note that existing models of wireless networks 
typically correspond to the case for which a(f) = 1 (or a positive constant independent of n) and a > 2q 
A common empirical model gives a(f) in dB/km for / in kHz as [16], [28]: 

P P 

l0loga(f) = a + a 1 p + a 2 -r—— FI + a 3 ———, (2) 

01 + J 2 b 2 + f 2 

3 A dense network [1], [6], [10] of unit area can also be considered as another fundamental network model, which will not be shown in 
this work. We remark that there exists either a bandwidth or power limitation (or both) according to the path-loss attenuation regimes in 
dense networks. 

4 The counterpart of a in wireless radio channels is the path-loss exponent. 



SUBMITTED TO IEEE TRANSACTIONS ON INFORMATION THEORY 



where {a , • ■ • ,03,61,62} are some positive constants independent of n. As mentioned earlier, we will 
allow the carrier frequency / to scale with n. Especially, we consider the case where the frequency scales 
at arbitrarily increasing rates relative to n, which enables us to really capture the dependence on the 
frequency in performance Jj The absorption a(f) is then an increasing function of / such that 

a(f) = 9 (e^ 2 ) (3) 



with respect to / for some constant c\ > independent of n. 

The noise rii at node i 6 {1, • - • , n} in an acoustic channel can be modeled through four basic sources: 
turbulence, shipping, waves, and thermal noise [16]. We assume that rii is the circularly symmetric complex 
additive colored Gaussian noise with zero mean and power spectral density (psd) N(f), and thus the noise 
is frequency-dependent. The overall psd of four sources decays linearly on the logarithmic scale in the 
frequency region 100 Hz - 100 kHz, which is the operating regime used by the majority of acoustic 
systems, and thus is approximately given by [16], [29] 

logiV(/)=a4-a5log/ (4) 

for some positive constants a 4 and a 5 independent of nu This means that N(f) = 0(1) since 

in terms of / increasing with n. From © and ©, we may then have the following relationship between 
the absorption a(f) and the noise psd N(f): 

N ^ - e (&[£»*) ■ <6) 

From the narrow-band assumption, the received signal y k at node k E {1, • • • , n} at a given time 
instance is given by 

Uk = ^ hkjXj + n k , (7) 

iei 



where 



(dJSm 



<k': 



V A (r ki ,f) 



(8) 



represents the complex channel between nodes i and k, Xi £ C is the signal transmitted by node i, and 
I C {1, • • • , n} is the set of simultaneously transmitting nodes. The random phases e> eki are uniformly 
distributed over [0, 2ti) and independent for different i, k, and time. We thus assume a narrow-band 
time-varying channel, whose gain changes to a new independent value for every symbol. Note that this 
random phase model is based on a far-field assumption [10], [11], [27]Jj which is valid if the wavelength 
is sufficiently smaller than the minimum node separation. 

Based on the above channel characteristics, operating regimes of the network are identified according 
to the following physical parameters: the absorption a(f) and the noise psd N(f) which are exploited 
here by choosing the frequency / based on the number n of nodes. In other words, if the relationship 
between / and n is specified, then a(f) and N(f) can be given by a certain scaling function of n from 
© and ©, respectively. 



'Otherwise, the attenuation parameter a(f) scales as 0(1) from {2j, which is not a matter of interest in this work. 

6 Note that in our operating frequencies, as = 1.8 is commonly used for the above approximation [16]. 

7 In [30], instead of simply taking the far- field assumption, the physical limit of wireless radio networks has been studied under certain 
conditions on scattering elements. Further investigation is also required to see whether this assumption is valid for underwater networks of 
unit node density in the limit of large number n of nodes. 



SUBMITTED TO IEEE TRANSACTIONS ON INFORMATION THEORY 5 

III. Cut-set Upper Bound 

To access the fundamental limits of an underwater network, a cut-set upper bound on the total throughput 
scaling is analyzed from an information-theoretic perspective [31]. Specifically, an upper bound based on 
the power transfer argument [10], [27] is established for extended networks. Note, however, that the 
present problem is not equivalent to the conventional extended network framework [10] due to different 
channel characteristics. Our interest is particularly in the operating regimes for which the upper bound is 
tight. 

Consider a given cut L dividing the network area into two equal halves as in [10], [27] (see Fig. [Q). 
Let Sl and D L denote the sets of sources and destinations, respectively, for the cut L in the network. 
More precisely, under L, source nodes Sl are on the left, while all nodes on the right are destinations 
Dl- In this case, we have an 0(n) x 9(n) multiple-input multiple-output (MIMO) channel between the 
two sets of nodes separated by the cut. 

In an extended network, we take into account an approach based on the amount of power transferred 
across the network according to different operating regimes, i.e., path-loss attenuation regimes. As pointed 
out in [10], [27], the information transfer from Sl to D L is highly power- limited since all the nodes in 
the set Dl are ill-connected to the left-half network in terms of power. This implies that the information 
transfer is bounded by the total received power transfer, rather than the cardinality of the set Dl- For 
the cut L, the total throughput T{n) for sources on the left is bounded by the (ergodic) capacity of the 
MIMO channel between Sl and D L under time-varying channel assumption, and thus is given by 



Tin) < maxE 

Ql>0 



logdet (l e{n) + ^L-H L Q L Hf ) 



(9) 



where 



where H L is the matrix with entries [H^]^ for i G Sl, k G D L , and Q L G C 6 ( n ) xQ ( n ) is the positive 
semi-definite input signal covariance matrix whose £>th diagonal element satisfies [Qzjfcfc < P for k G Sl- 
The relationship © will be further specified in Theorem [Q Before that, we first apply the techniques 
of [10], [32] to obtain the total power transfer of the set Dl- These techniques involve the relaxation of 
the individual power constraints to a total weighted power constraint, where the weight assigned to each 
source corresponds to the total received power on the other side of the cut. To be more specific, each 
column i of the matrix H L is normalized by the square root of the total received power on the other side 
of the cut from source i G Sl. From (OQ) and ([8]), the total power PJ?' received from the signal sent by 
the source i is given by 

if = M 4 \ do) 

C °ken L 

for some constant c > independent of n. For convenience, we now index the node positions such 
that the source and destination nodes under the cut L are located at positions (—i x + 1, i y ) and (k x , k y ), 
respectively, for i x , k x — 1, • • • , y/n/2 and i y , k y = 1, • • • , \fn. The scaling result of dl defined in (fTT|) 
can then be derived as follows. 

Lemma 1: In an extended network, the term d/£ in (fTTI) is 

df = 6 {i l - a a(f)-^) , (12) 

where — i x + 1 represents the horizontal coordinate of node i G Sl for i x = 1, • • • , y/n/2. 
The proof of this lemma is presented in Appendix [A] The expression © is then rewritten as 



maxE 

Qi>o 



logdet (le (n ) + j^FlQlfA 



(13) 
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where F L is the matrix with entries [F L ]ki 
Dl- Here, Ql is the matrix satisfying 

Ql 



= [Hx,]fci, which are obtained from (fTTj) . for i G Sl, k G 



- ki 



,(k)S) 



dfdl> [Q L 



ki - 



,(') 



which means that tr(Q L ) < Xlies ^l' ( et l ua l to the sum of the total power received from each source). 

We next examine the behavior of the largest singular value for the normalized channel matrix F L , 
and then show how much it affects an upper bound on (fl"3l) . We first address the case where F^ is 
well-conditioned according to the attenuation parameter a(f). 

Lemma 2: Let F L denote the normalized channel matrix defined by the expression (fl"3T) . Under the 
attenuation regimes a(f) = fi ((1 + e )^) for an arbitrarily small e > 0, we have that 

I 21 



E[\\F L \\l] <c 2 log 



n 



(14) 



for some constant c 2 > independent of n. 

The proof of this lemma is presented in Appendix [Bj Note that the matrix F^ is well-conditioned as 
a(f) scales exponentially with respect to y/n (or faster). Otherwise, i.e., if a(f) = o ((1 + eo)^), the 
largest singular value of F L scales as a polynomial factor of n, thus resulting in a loose upper bound on 
the total throughput. Using Lemma |2l we obtain the following result. 

Lemma 3: Under a(f) = fi ((1 + eo)^), the term (fT3l) is upper-bounded by 



n 



mr> 



E^ 



(i) 

L 



(15) 



i&S L 



>« 



for arbitrarily small positive constants e and e, where P L is given by (flOl) . 
Proof: Equation (TT3l) is bounded by 



m&xE 

Ql>0 

+ m&xE 

Ql>0 



log det I I e („) + 
log det ( I e („) + 



H 



N(f) 
1 



FlQlF l 
F L Q L Ff 



l £F r 



(16) 



Here, the event £f l refers to the case where the channel matrix F^ is accidently ill-conditioned and is 
given by 

£ FL = {\\F L \\l>n<} 

for an arbitrarily small constant e > 0. Suppose that a(f) — Q ((1 + eo)^). Then, by using the result of 
Lemma [2] and applying the proof technique similar to that in Section V of [10], it is possible to prove 
that the first term in (fT6l) decays polynomially to zero with arbitrary exponent as n tends to infinity, and 
for the second term in (fT6l) . it follows that 



max E 

Ql>0 



log det ( I 6 („) 
1 



N{f) 



F L Q L Ff )l 



< ma.xE 

Ql>0 

c 2 log n 



N(f) 



\F L \\' 2 tTlQ L )l £ . 



< 



< 



maxtr Q L 
^{J) Ql>0 



n 



N(S) ■ 



Up. 



(0 

L 



i&S L 
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for some constant c 2 > independent of n, where the second inequality holds by Lemma [2] This completes 
the proof of this lemma. ■ 

Note that (fl"5l) represents the total amount of received signal-to-noise ratio from the set Sl of sources 
to the set Dl of destinations for a given cut L. We are now ready to show the cut-set upper bound in 
extended networks. 

Theorem 1: For an underwater regular network of unit node density, where the absorption coefficient 
a(f) scales as O ((1 + e )^) for an arbitrarily small e > 0, the total throughput T{n) is upper-bounded 
by 

T{n) * irmn <17) 

where c 3 > is some constant independent of n and e > is an arbitrarily small constant. 

Proof: Suppose that a(f) — Q ((1 + eo)^). Then from Lemmas Q] and [3] we obtain the following 
upper bound on the total throughput T(n): 



Tin) < — -r-J^Pd 



L 



Pn e 



Vn/2 Jh 



— EE4" 



< c t Pn l ' 2 +' ^" 



'L 



n /) .4; ir l <f) 



< Cl PnV^ ■ " : 



y — 



< 






a(/)^(/)' 

where c 4 and c 5 are some positive constants independent of n, which is equal to (fT7l) . This completes the 
proof of the theorem. ■ 

Note that this upper bound is expressed as a function of the absorption a(f) and the noise psd N(f) 
while an upper bound for wireless radio networks depends only on the constant value a [10]. We remark 
that when a(f) = o ((1 + eo)^"), the upper bound becomes boosted by a certain polynomial factor of n 
(up to 0(\/n)) compared to the case shown in (fi7])l 8 l In addition, using © in (fTTT) results in 

r(n) = o l(loga(/)r5/2nl/2+£ 



for some constant a 5 > shown in ©. Finally, another expression for the condition in which the upper 
bound in (TT7T ) holds is shown as follows. 

Remark 1: We examine the relationship between the carrier frequency / and the number n of nodes 
such that our upper bound holds. By using ([3]) and the regimes a(f) = fi((l + e )v™), we can also obtain 
the following condition: 

which means that if / scales faster than n l / A , then the result in (IT7T) is satisfied. 

8 This statement could be rigorously proved by following the steps similar to those shown in Lemmas [2] and [3] even if the details are not 
shown in this paper. 
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IV. ACHIEVABILITY RESULT 

In this section, we show that the considered transmission scheme, commonly used in wireless radio 
networks, is order-optimal in underwater networks. Under a regular network of unit node density, the 
conventional MH transmission is used and its achievable throughput scaling is analyzed to show its order 
optimality. 

The nearest-neighbor MH routing protocol [1] will be briefly described with a slight modification. The 
basic procedure of the MH protocol under our extended regular network is as follows: 

• Divide the network into square routing cells, each of which has unit area. 

• Draw a line connecting a S-D pair. A source transmits a packet to its destination using the nodes in 
the adjacent cells passing through the line. 

• The full power is used, i.e., the transmit powers at each node is P. 

Instead of original (continuous) MH transmissions, a bursty transmission scheme [10], [27], which uses 
only a fraction l/a(f)N(f) of the time for actual transmission with instantaneous power a(f)N(f)P 
per node, is used to simply apply the analysis for networks with no power limitation to our network 
model. With this scheme, the received signal power from the desired transmitter, the noise psd, and the 
total interference power from the set I C {1, • • • ,n} have the same scaling, i.e., Q(N(f)), and the 
(instantaneous) received signal-to-interference-and-noise ratio (SINR) is kept at 0(1) under the narrow- 
band model, which will be obviously shown in the proof of Theorem [2] 

The achievable rate of MH is now shown by quantifying the amount of interference. 

Lemma 4: Suppose that a regular network of unit node density uses the MH protocol with burstiness. 
Then, the total interference power from other simultaneously transmitting nodes, corresponding to the 
set / C {1, ■ ■ • , n}, is upper-bounded by 0(iV(/)), where N(f) denotes the psd of noise n, at receiver 
i 6 {l,- ,n}. 

Proof: There are 8k interfering routing cells, each of which includes one node, in the A;-th layer 4 
of the network as illustrated in Fig. [2] Then from (OQ), ©, and ([8]), the total interference power at each 
node from simultaneously transmitting nodes is upper-bounded by 



Y^ ( h Mf)N(f)P 8N(f)P ^ 1 



k=i 



oo 
2-^ „( f\k-l 



< BN(f)P 
~ co ^aifY 

< ceN(f), 

where c and cq are some positive constants independent of n, which completes the proof. ■ 

Note that the signal power no longer decays polynomially but rather exponentially with propagation 
distance in our network. This implies that the absorption term a(f) in CQ) will play an important role in 
determining the performance. It is also seen that the upper bound on the total interference power does 
not depend on the spreading factor a. Using Lemma H] it is now possible to simply obtain a lower bound 
on the capacity scaling in the network, and hence the following result presents the achievable rates under 
the MH protocol. 

Theorem 2: In an underwater regular network of unit node density, 

T(n) = a \AfWU)) <18) 

is achievable. 

Proof: Suppose that only a fraction l/a(f)N(f) of the time for actual transmission is used under 
the MH protocol with burstiness. Then, the SINR seen by any receiver is expressed as f2(l) with an 
instantaneous transmit power a(f)N(f)P since the total interference power is given by 0(N(f)). Since 
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the Gaussian is the worst additive noise [33], [34], assuming it lower-bounds the throughput. Hence, by 
assuming full CSI at the receiver, the achievable throughput per S-D pair is lower-bounded by 

W M1 + SINR) 

1 . A N(f)P/co \ 



~ a(f)N(f) l ° g V/ ' N(f) + c 6 N(f)J 
for some positive constants cq and cq independent of n, thereby providing the rate of 

1 \ 



n 



a(f)N(f)J- 



Since the number of hops per S-D pair is given by 0{y/n), there exist £l(y/n) source nodes that can 
be active simultaneously, and therefore the total throughput is finally given by (fl"8l) . which completes the 
proof of the theorem. ■ 

Now it is examined how the upper bound shown in Section [Till is close to the achievable throughput 
scaling. 

Remark 2: Based on Theorems [Q and [2l when a(f) = fi ((1 + eo)^), i.e., / = Vtin 1 ^), it is easy 
to see that the achievable rate and the upper bound are of the same order up to n e , where e and eo are 
vanishingly small positive constants. MH is therefore order-optimal in regular networks with unit node 
density under the above attenuation regimes. 

We also remark that applying the hierarchical cooperation strategy [10] may not be helpful to improve 
the achievable throughput due to long-range MIMO transmissions, which severely degrade performance 
in highly power-limited networks^ To be specific, at the top level of the hierarchy, the transmissions 
between two clusters having distance 0{y/n) become a bottleneck, and thus cause a significant throughput 
degradation. It is further seen that even with the random phase model, which may enable us to obtain 
enough degrees-of-freedom gain, the benefit of randomness cannot be exploited because of the power 
limitation. 

V. Extension to Random Networks 

In this section, we would like to mention a random network configuration, where n S-D pairs are 
uniformly and independently distributed on a square. 

We first discuss an upper bound for extended networks of unit node density. A precise upper bound can 
be obtained using the binning argument of [10] (refer to Appendix V in [10] for the details). Consider 
the same cut L, which divides the network area into two halves, as that in the regular network case. For 
analytical convenience, we can artificially assume the empty zone El, in which there are no nodes in the 
network, consisting of a rectangular slab of width < c < / = 1 1/4 , independent of n, immediately to the 

right of the centerline (cut), as done in [27] (see Fig. |3]o Let us state the following lemma. 

Lemma 5: Assume a two dimensional extended network where n nodes are uniformly distributed. When 
the network area is divided into n squares of unit area, there are fewer than log n nodes in each square 
with high probability. 

Since the result in Lemma [5J depends on the node distribution but not the channel characteristics, the 
proof essentially follows that presented in [4]. By Lemma [5J we now take into account the network 
transformation resulting in a regular network with at most log n and 2 log n nodes, on the left and right, 

9 In wireless radio networks of unit node density, the hierarchical cooperation provides a near-optimal throughput scaling for the operating 
regimes 2 < a < 3, where a denotes the path-loss exponent that is greater than 2 [10]. Note that the analysis in [10] is valid under the 
assumption that a is kept at the same value on all levels of hierarchy. 

10 Although this assumption does not hold in our random configuration, it is shown in [27] that there exists a vertical cut such that there 
are no nodes located closer than < c < / - 1 1/4 on both sides of this cut when we allow a cut that is not necessarily linear. Such an 
existence is proved by using percolation theory [4], [35]. This result can be directly applied to our network model since it only relies on 
the node distribution but not the channel characteristics. Hence, removing the assumption does not cause any change in performance. 
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respectively, at each square vertex except for the empty zone (see Fig. [3]). Then, the nodes in each square 
are moved together onto one vertex of the corresponding square. More specifically, under the cut L, the 
node displacement is performed in the sense of decreasing the Euclidean distance between source node 
i E Sl and the corresponding destination k G Dl, as shown in Fig. [3l which will provide an upper 
bound on dl in (fTTT) . It is obviously seen that the amount of power transfer under the transformed regular 
network is greater than that under another regular network with at most log n nodes at each vertex, located 
at integer lattice positions in a square region of area n. Hence, the upper bound for random networks is 
boosted by at least a logarithmic factor of n compared to that of regular networks discussed in Section Hill 
Now we turn our attention to showing an achievable throughput for extended random networks. In 
this case, the nearest-neighbor MH protocol [1] can also be utilized since our network is highly power- 
limited. Then, the area of each routing cell needs to scale with 2 log n to guarantee at least one node in a 
cell [1], [6]lH| Each routing cell operates based on 9-time division multiple access to avoid causing large 
interference to its neighboring cells [1], [6]. For the routing with continuous MH transmissions (i.e., no 
burstiness), since per-hop distance is given by 9 ( A/log n), the received signal power from the intended 
transmitter and the SINR seen by any receiver are expressed as 

c 7 P 

(logn) a / 2 a(/) C8 v / ^g^ 

and 

n 



(logn) a / 2 a(/) c 8v / ^JV(/), 

respectively, for some constants c 7 > and eg > \/2 independent of n. Since the number of hops per 
S-D pair is given by 0(^n/\ogn), there exist Q,(\Jn/ 'logn) simultaneously active sources, and thus the 
total achievable throughput T(n) is finally given by 

„l/2 

Tin) = Q ' 



(log nY a+1 )/ 2 a(f)^^^N{f) , 

for some constant cs > v2 independent of n (note that this relies on the fact that log(l + x) can be 
approximated by x for small x > 0). Hence, using the MH protocol results in at least a polynomial 
decrease in the throughput compared to the regular network case shown in Section ITvPl This comes from 
the fact that the received signal power tends to be mainly limited due to exponential attenuation with 
transmission distance Q(y/\ogn). Note that in underwater networks, randomness on the node distribution 
causes a huge performance degradation on the throughput scaling. Therefore, we may conclude that the 
existing MH scheme does not satisfy the order optimality under extended random networks regardless of 
the attenuation parameter a(f). 

VI. Conclusion 

The attenuation parameter and the capacity scaling laws have been characterized in a narrow-band 
channel of underwater acoustic networks. Provided that the carrier frequency / scales at arbitrary rates 
relative to the number n of nodes, the information-theoretic upper bound and the achievable throughput 
were derived as a function of the attenuation parameter a(f) in extended regular networks. Specifically, 
based on the power transfer argument, the upper bound was shown to decrease in inverse proportion 
to a(f). In addition, to show the achievability result, the nearest-neighbor MH protocol was introduced 

"When methods from percolation theory are applied to our random network [4], [35], the routing area constructed during the highway 
phase is a certain positive constant that is less than 1 and independent of n. The distance in the draining and delivery phases, corresponding to 
the first and last hops of a packet transmission, is nevertheless given by some constant times log n, thereby limiting performance, especially 
for the condition a(f) — w(l). Hence, using the protocol in [4] indeed does not perform better than the conventional MH case [1] in random 
networks. 

12 In terrestrial radio channels, there is a logarithmic gap in the achievable scaling laws between regular and random networks [1], [25]. 
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with a simple modification, and its throughput scaling was analyzed. We proved that the MH protocol is 
order-optimal as long as the frequency / scales faster than n 1//4 . Our scaling results were also extended 
to the random network scenario, where it was shown that the conventional MH scheme does not satisfy 
the order optimality for all the operating regimes. 

Suggestions for further research include, in dense networks of unit area, analyzing an upper bound and 
designing an achievable scheme whose throughput scaling is close to the upper bound. 

Appendix 

A. Proof of Lemma \B 

Upper and lower bounds on d£ can be found by using the node-indexing and layering techniques similar 
to those shown in Section VI of [32]. As illustrated in Fig. HI layers are introduced, where the ?-th layer l[ 
of the network represents the ring with width 1 drawn based on a source node, whose coordinate is given 
by (—i x + 1, i y ), where i E {1, • - • , \/n}. More specifically, the ring is enclosed by the circumferences of 
two circles, each of which has radius i x + i — l and i x + i — 2, respectively, at its same center (see Fig. HI). 
We can see that there exist Q(i x + i) nodes in the layer l[ since the area of /• is given by n(2i x + 1i — 3). 
Then from (fTTI) . the term d£ is given by 



'x 

yfn/2 -y/n 

) _ l 
b L 



^=ytt 



c ° f i £i «i, + fe - l) 2 + («, - W' 2 a(/)V(-+t.-i) 2 +('.-*.) 2 



It is further assumed that all the nodes in each layer are moved onto the innermost boundary of the 
corresponding ring, which provides an upper bound for d[ . Since there is no node, located on the right 
half of the cut L, in the first layer l[, d^ is then upper-bounded by 



4"< 



1 oo 



Cgk' 



< 



Co fc tt k ' aa W 

oo .. 

eg v^ 1 



wr 1 fe tt a{f)k ' 

< ^^ (_l_ r i d 



< 



ClO 



where Co, eg, and cio are some positive constants independent of n. Here, the fourth inequality holds 
since a(f) > 1. To get a lower bound for d[ , nodes in each layer are now moved onto the outermost 
boundary of the corresponding ring. Let 5i denote the fraction of nodes that are placed on the right half 
of the network among nodes in the i-th layer l[, which is obviously independent of n irrespective of 
i E {1, • • • , y/n/2}. Thus, the lower bound similarly follows 

^+vV2-l 



''L 



i/) . -i y" c 9h'-i x +ik' 

i 

i x +y^n/2-l 



> 



c ^ k ,a a(f) k ' 

U k'=i x W ! 

c 9 min{5i, ••• ,5^/2} 



fc tt k ' a -^(f) k ' 



> 



Cll 



where c , c 9 , and en are some positive constants independent of n, which finally yields (fT2l . This completes 
the proof. 



SUBMITTED TO IEEE TRANSACTIONS ON INFORMATION THEORY 



12 



B. Proof of Lemma [2] 

Since the size of matrix F^ is given by 0(n) x Q(n), the analysis essentially follows the argument 
in [10] with a slight modification (refer to Appendix III in [10] for more precise description). Suppose 
that a(f) = fi ((1 + e )"^) for an arbitrarily small e > 0. Then in the following, from the result of 
Lemma [Q we derive J2keD L \i F L}ki\ 2 and an upper bound for J2iss L \[ F L]ki\ 2 - 

2 



keD L 



k£D L 



fH 



L \ki 



E 



k&D L 



fH 



L\ki\ 



zZk&Dr. A( r ki, f) 



1. 



where the second equality comes from ©, ([8]), and (fill) , and 

2 



ies L 



1 



fH 



ijfci 



= E 

ies L y^ L 

tJtL/2 y/n 

<c yy 

£ri £ri ((4 + fc x - 1)2 + (^ - ^) 2 ) Q/2 a(/)v / ^+ fc - 1 ) 2 +fe-^) 2 

< c °EE 



ir 1 a(/) 1 



=1 h,=1 \/(^ + kx~ I) 2 + (« 
v^/2 1 / VH 

_ Z_T + fc^ 11. 



o(/)*- 



1 



V^/2 , , 






^i o(/)V«+?-^ 

1 



< Co 



£=i »x I ^ a ( / )^/(v / ^f+^) 



v^/2 1 / v 7 ^ -j 

^ ^ I ^ (1 + ec 



1C12«S 



ix=l 



< Co I 1 + 

< ci 3 log n 



v^/2 l 

—dx 



x 



eoj 



+ 



\/n 



:i+e ; 



ci2y 



dy 



for some positive constants Cyi and C13 independent of n, where the second inequality holds since a = 1 
provides the highest upper bound for all 1 < a < 2. The fifth inequality comes from the fact that 
a(f) = Cl ((1 + e )v"). Hence, it is proved that both scaling results are the same as the regular network 
case shown in [10]. 

Now we are ready to prove the inequality in (fT4l) . Following the same line as that in Appendix III 
of [10], we thus have 



£[tr((F£F L ) 9 )] <C>(c 14 log 



n 



where C a 



(2?)! 



tM ; ^, is the Catalan number for any q and c i4 > is some constant independent of n. Then, 
from the property ||F L ||| = lim^oo tr((F^F L )' 7 ) 1 /' ? (refer to [36]), the expectation of the term ||F L ||| is 
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upper-bounded by 



E[]\F L ]\l] <ta(B[tr((Ff Fl )')]) 

< lim {C q n (C14 log n 

q— >oo 

= 4ci 4 logn, 






where the equality holds since lim^oo C q = 4. Here, the first inequality comes from dominated 
convergence theorem and Jensen's inequality. This completes the proof. 
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Fig. 1. The cut L in a two-dimensional extended regular network. Si, and Dl represent the sets of source and destination nodes, respectively. 
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Fig. 2. Grouping of interference routing cells in extended networks. The first layer l\ represents the outer 8 shaded cells. 
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Fig. 3. The node displacement to square vertices, indicated by arrows. The empty zone El with width constant c is assumed for simplicity. 
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Fig. 4. Grouping of destination nodes in extended networks. There exist Q(i x ) nodes in the first layer l[. This figure indicates the case 
where one source is located at the position (— i x + l,i y ). The destination nodes are regularly placed at spacing 1 on the right half of the 
cut L. 



